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We calculate one-loop matching factors for bilinear operators composed of improved staggered 
fermions. We compare the results for different improvement schemes used in the recent literature, 
all of which involve the use of smeared links. These schemes aim to reduce, though not completely 
'^v^' ' eliminate, O(a^) discretization errors. We find that all these improvement schemes substantially 

, reduce the size of matching factors compared to unimproved staggered fermions. The resulting 

■ corrections are comparable to, or smaller than, those found with Wilson and domain-wall fermions. 

^N) ' In the best case ("Fat- 7" and mean-field improved HYP links) the corrections are 10% or smaller 

! at 1/a = 2 GeV. 

o . 

O: 

, . I. INTRODUCTION 

(N : 

Improved staggered fermions are an attractive choice for numerical simulations of unqucnched QCD [1-11]. They 
CO ' maintain the positive features of unimproved staggered fermions — smaller CPU requirements than other fermion dis- 
^ ■ cretizations, a remnant chiral symmetry, and 0{a?) discretization errors — while potentially avoiding its drawbacks — 
large flavor symmetry breaking and large perturbative matching factors. We have begun a program of calculations 
. of electroweak matrix elements using these fermions, and thus need to decide between the different improvement 
QQ ' schemes that have been suggested in the recent literature. Since we intend at first to use one-loop perturbation theory 
, to match lattice and continuum operators, it is important that the matching factors for relevant operators are close 
C^l ' to unity. In this paper we calculate these matching factors for the bilinear operators which form the building blocks 
\ of the four-fermion operators we intend to use in our matrix element calculations. We expect our results to be a 
good guide to the size of corrections for the four-fermion operators themselves — typically the one-loop contributions 
get roughly doubled. In any case, finding small one-loop corrections for bilinears is a prerequisite for proceeding to 
JL ■ four-fermion operators. 

We stress that we are using the term "improvement" loosely in this work. Although the improved actions and 
, operators that we use are motivated by the Symanzik program, we are not following this program systematically. This 
would involve improving the gauge and fermion actions, and the operators, so as to remove O(a^) errors either order 
by order in perturbation theory or non-perturbatively. Such a program is much more difficult for staggered fermions 
' than for Wilson-like fermions, particularly with the operators we use which are spread out over a 2* hypercube.^ Our 
?H ] aim is to change the fermion action and operators such that the 0{a?) corrections are reduced from the large size 
typical of unimproved staggered fermions to the size seen with Wilson, Domain wall or overlap fermions. We do not 
improve the gauge action. 

The plan of this paper is as follows. In the next section we describe the alternatives we have considered for improved 
operators and actions. Section III collects the new features of the Feynman rules that are introduced by improvement. 
In sec. IV we present analytic results for the one-loop matching constants. We then, in sec. V, describe how to do a 
second level of mean-field improvement of the bilinear operators. We close in sec. VI with the numerical results and 
a discussion of their implications. We collect some definitions in an Appendix, along with the results that allow us to 
push the analytic calculation one step further than in previous work. 

We lean heavily on the notation and methodology of Ref. [14], and for brevity we refer to that paper as PS in the 
following. 



'^For recent work on further improving staggered fermions see Refs. [12,13]. 
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II. IMPROVED ACTIONS AND OPERATORS 



Lepage [2,8] and Lagae and Sinclair [6] have argued that flavor symmetry breaking can be substantially reduced 
by suppressing the coupling of high momentum gluons which connect "physical" quarks residing at different corners 
of the Brillouin zone. Such suppression is also expected to reduce the size of one-loop contributions to perturbative 
matching factors, since, as noted by Golterman [15], their size is largely due to tadpole-type diagrams which involve 
these flavor-changing vertices. 

The flavor-changing coupling can be suppressed by replacing the standard "thin" link with some form of "smeared" 
link in the quark covariant derivative. Various options have been tried, and we consider here two choices which have 
been successful at reducing flavor symmetry breaking in pion masses, are relatively local, and have been extensively 
studied: the "Fat- 7" link introduced by Orginos and Toussaint [7] and studied numerically by the MILC collabora- 
tion [10], and the HYP link introduced by Hascnfratz and Knechtli [16]. We refer to the original papers for the details 
of the constructions and do not repeat them here. Both are "fattened" by averaging over paths containing links in 
some or all of the transverse directions (and which in the Fat-7 case are up to 7 links long), and in this way they 
reduce the coupling to gluons with transverse momenta of 0{l/a). In some sense the Fat-7 link is the simplest choice 
which accomplishes this, while the HYP link involves an average over more paths. On the other hand, the HYP link 
involves three levels of APE-like smearing with projection back into SU(3). Simulations show that such smearing 
is very effective at reducing flavor-symmetry breaking. Indeed, using the HYP links leads to a greater reduction in 
flavor symmetry breaking than the Fat-7 links. 

The introduction of smeared links can be viewed as one part of the Symanzik improvement program applied 
at tree-level to staggered fermions. Complete removal of O(a^) terms from fermion vertices requires two other 
improvements [8]. First, the smearing of the links introduces an O(a^) correction to the flavor-conserving quark-gluon 
coupling. This can be removed by a adding to the smearing a 5-link "double staple" - we refer to this as the "Lepage 
term". Second, the O(a^) corrections to the fermion propagator need to be removed, and this can be done by adding 
a next-to-nearest-neighbor derivative, the "Naik term" [19]. 

An observation of practical relevance is that the Naik term is the only part of the improvement of the fermion 
vertex that cannot be accomplished simply by changing the links in the unimproved staggered action. In other words, 
if one does not include the Naik term, and if one is interested in calculating propagators on configurations that have 
been already been generated (whether quenched or unquenched) then the practical implementation of smeared links 
is simple: one calculates the smeared links, and then uses an unimproved staggered inverter. 

Complete tree-level O(a^) improvement of physical quantities requires, in addition to the improvement of fermion 
vertex outlined above, the use of a tree-level (or more highly) improved gauge action. The previous discussion implies, 
however, that 0{a^) errors from the gauge action are not responsible for the large flavor-symmetry breaking or the 
large "tadpole" contributions to matching factors. 

With these general comments in mind, we can now explain our choices of action and operators. We use the single- 
plaquette Wilson gauge action, since this is the action we are using in our present simulations.^ For the fermion 
action, we keep the original staggered form (without the Naik term), but use various types of smeared links. The 
only exception is case (5) below, in which we keep the Naik term. Finally, for the bilinears we use the standard 
hypercube form (the definition of which is given in the Appendix), rendered gauge invariant by including the average 
of the product of links along the shortest paths between the quark and anti-quark field. ^ The only improvement of 
the operators that we consider is the use of smeared links. Intuitively, the reduction in fluctuations in these links will 
reduce the flavor symmetry breaking between bilinears [16] . In all cases we use the same type of smeared links in the 
operators as in the action, so that the hypercube vector current is conserved [except in case (5)]. 

The specific choices of links we consider are: 

1. The original gauge links, tadpole improved (following the prescription of Rcf. [18] as implemented in PS). We 
use the fourth root of the average plaquette to determine the "average link" mq. This yields (tadpole- improved) 
unimproved staggered fermions and unimproved operators, and allows us to check our results against those in 
PS. 



■^For this reason, we cannot compare our results with those of Rcf. [20], since these authors use an improved gauge action. 

^We do not consider here so-called "Landau gauge" operators — those rendered gauge invariant by transforming to Landau 
gauge and then leaving out the links. These are not useful for matrix elements involving "eye" diagrams, because they allow 
mixing with lower-dimension gauge non- invariant operators [17]. They are also subject to uncertainties due to the presence of 
Gribov copies. 
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2. Fat-7 smeared links, built out of tadpole improved links (as in the numerical implementation of Ref. [7,10]). We 
stress again that we use these smeared links both in the action and in the bilinear operators. 

3. Fully O(a^) improved smeared links, i.e. Fat-7 links with the Lepage "double-staple" term added, again both 
in the action and the operators. 

4. Links smeared according the HYP prescription of Hasenfratz and Knechtli [16] , again both in the action and 

the operators. Three parameters, Q!i_3, need to be specified to completely define HYP smearing, and we focus 
on two choices, as described below. We also consider a variant in which we tadpole improve the smeared links 
themselves (section V). 

In addition, we consider a final choice of action and operators: 

5. Following the "Asqtad" action introduced by Lepage [8] and used extensively by the MILC collaboration [10], 
we add the Naik term to the action of choice 3, while taking the same operators as in choice 3. In the Naik 
term alone, we use the original unsmeared gauge links (tadpole improved). 

This differs from the "Asqtad" action, however, because we use the unimproved Wilson gauge action, whereas "Asqtad" 
includes an improved gauge action. We thus refer to our choice as the "Asqtad-like" action. Our expectation is that 
the choice of gauge action has relatively little impact on the size of matching factors, and particularly on the variation 
of these factors between bilinears having the same spin and different flavor. 

III. FEYNMAN RULES 

The Feynman rules for unimproved staggered fermions are standard. In the notation we use here, they can be found 
in PS, and we do not repeat them. We discuss only the changes introduced by smearing the links and including the 

Naik term. 

We consider first the effect of smearing the links. For all except tadpole diagrams (i.e. those in which two gluons 
emerge from a single, possibly smeared, link), the only effect is to change the coupling to the underlying gluon field. 
With the unimproved action, a link in the ^-th direction couples only to a gluon A^{k) with v = ^. The smeared 
links, however, couple to Aii(k) for all v, and the extra factor this introduces can be conveniently written as 

5u,t.D^{k) + {l~5,,^)GuAk)- (1) 
The diagonal and off-diagonal couplings can be decomposed, respectively, as 

D^{k) = l-diJ2sl + d2 J2 - d^slspl - d4 (2) 
with Sv = sin(fcy/2), etc., and 

a,^(fc) =s^s.G.,^(fc) (3) 

G^Ak) =di-d2 2 + + d^sl , (4) 

where all indices (/x, p, a) are different. 

The coefficients di_4 distinguish the different choices of links: 

1. Unimproved (choice 1 above): 

di =0, d2 = 0, d3 = 0, di = 0. (5) 

2. Fat-7 links (choice 2 above): 

di = 1, d2 = 1, ^3 = 1, d4 = 0. (6) 

3. O(a^) improved links (choices 3 and 5 above): 

=0, d2 = 1, ds = 1, di = 1. (7) 
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4. HYP smeared links (choice 4 above): 



rfi = (2/3)ai(l + a2(l + Q!3)), ^2 = (4/3)aia2(l + Zag), da = 8aia2a3, ^4 = 0. (8) 

We consider two choices for the ai. The first was determined in Ref. [16] using a non-perturbative optimization 
procedure: ai = 0.75, a2 = 0.6 as = 0.3. This gives 

c;i=0.89, ^2 = 0.96, ^3 = 1.08, ^4 = 0. (9) 

The second is chosen so to remove O(a^) flavor-symmetry breaking couplings at tree level. This choice, ai = 7/8, 
C(2 = 4/7 and as = 1/4, gives 

di = 1, d2 = 1, = 1, di = 0, (10) 

i.e. the same as for Fat-7 links. 

These results agree with those of Refs. [20,21], but are written here in a somewhat different notation. The fact that all 
four choices can be collected in this form simplifies the resulting one-loop calculations. It is particularly noteworthy 
that the Fat-7 and HYP vertices can be made identical, showing that these two choices cannot be distinguished by 
their flavor breaking effects in perturbation theory [21]. The one-loop matching factors for these two choices are not, 
however, identical, because the tadpole contributions differ. 

For tadpole diagrams, which involve two-gluon vertices, the differences between the actions are more complicated, 
and will be given explicitly below. 

The inclusion of the Naik term alters the Feynman rules in several ways. In the fermion propagator, all factors of 
Sfj, = sinp^ are replaced: 

S^^S^^ =S^{l + dNsl/&). (11) 

Here wc have introduced a fifth coefficient d{^ which distinguishes the different choices of action: d^ = unless the 
Naik term is included, in which case djv = 1- This device allows us to write most of our results in a way which holds 
for all choices of action and operators. 

The one-gluon vertex is also changed by the Naik term, but this can only be represented in a simple way if one 
of the quarks in the vertex has vanishing physical momentum [k^ = (0,7r/a)]. In this case, the diagonal part of the 
vertex changes as follows: 

D^^D'll = D^ + dNsl/&. (12) 
This substitution works for all except the self-energy diagram, which we consider explicitly below. 



IV. ANALYTIC RESULTS FOR MATCHING CONSTANTS 



The one-loop matching relations take the general form, 

0co„t ^ 0iat +Cp^ Y.^kj2di In(Ma) + Ci,)Of (13) 

j 

where Cp = 4/3 is the color Casimir factor, fi is the renormalization scale of the continuum operators, and i and j 
run over all the different possible bilinears in a four-flavor theory. The explicit forms of the operators are given in the 
Appendix. The constants di are proportional to the one-loop anomalous dimensions of the bilinears, 7^°^ = —2CFdi. 
They depend only on the spin of the bilinear, and are di = (3, 0, —1) for spins {S/P, V/A, T). The finite part of the 
coefficient can be written 

Cij = 6ij [rfj(7B - Foooo) + ts] - clj , (14) 

with ts depending on the continuum renormalization scheme. For the NDR scheme ts = (—0.5,0,1.5) for spins 
{S/P, V/A, T). The conversion to other schemes is given in PS. The constants are 7^ = 0.577216 and Fqooo = 4.36923. 
Finally, the "lattice" part of the coefficient can be broken up as follows: 

= Xij + % {Y, + T, + Z) . (15) 

Here X, Y, T and Z refer to contributions from the different types of diagrams using the notation of PS, as illustrated 
in Fig. 1. This equation incorporates the fact, derived below, that only the X diagrams lead to mixing among bilinears. 
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FIG. 1. Notation for diagrams contributing to matching factors. 

All our calculations are done in Fcynman gauge. We have checked our results by doing two independent calculations 
using different methods — the first following PS and the second using the methods presented in Refs. [22,23]. 

A. X diagrams 

The calculation follows the same steps as in PS, except for two changes. 

• We have been able to carry out the calculation analytically apart from one final integral, using the results 

(A8,A9) given in the Appendix. 

• The improved vertex eq. (1) allows propagation from a smeared link in the ji-ih. direction to another in any 
direction, even in Feynman gauge. It is useful to distinguish between the case where the second link is also in 
the /x-th direction, for which the gluon propagator is multiplied by 

P;^^{k)=D^^{kf + Y,G.Ak)\ (16) 

and the case where the second link is in a different direction p ^ /x, for which the multiplying factor is 

D^^ G^,, + D^Gp,^ + J2 G.,^G,,p ^ ^s^SpO^Jik) . (17) 

The superscripts emphasize the fact that there is a possible Naik term at both ends of the propagator. 
Using these results, we find the following expression for the diagonal part of the contributions from the X diagrams: 



E 



fc L 



(-1) 



/ v<..si)'o;rBFV,(fc) \i-{-if^+'' 



(18) 



Here = IBtt^ /-t ''M — cos(A;^/2) and = sm{kfj,). The functions arising from boson and fermion 
propagators are, respectively. 
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m = TT^, F{k) = ^^-^. (19) 



For the sake of brevity we do not show the argument k of these functions or of P^'^ and Oj^J^ in eq. (18) and in the 
following. The index i in eq. (18) labels the spin and flavor of the operator, and the /A term on the first line is the 
conventional integral used to cancel divergences. The function Vj(fc) is defined in eq. (AS); we stress again that the 
use of this equation leads to a simpler form than that given in PS. 

Only the X diagrams lead to mixing, i.e. non-zero values for X^ , i ^ j- We find that we can also give explicit 
expressions for the mixing terms using eq. (A9). As for unimproved staggered fermions. the mixing that occurs at 
one-loop turns out to be only a subset of that allowed by the hypercubic symmetry group. The non-zero mixing 
coefficients are (using the definitions in PS — see table II) 

CVVM = / ABF SiSi S2S2 (is C3 - S3) > (20) 

J k 

cvAM = f ABFHrs'^S2S^{-P^''clc3Ci + Ofa'^s^pca - C3]c4) , (21) 

Jk 

CVTM = f 2BFh,S^S2S^{-Pi'''cl[C3 + C4] + O^/" sl[2c2 -CS + c4) , (22) 
J k 

CTAM = [ 2BFhiS^S2S^{Ps''''cl[c^ - Cs] + 0^2''SI[2C2 - C3 - C4]) , (23) 
J k 

where = cos(A;^). 

B. Y diagrams 

Y diagrams involve the gluon connecting an external quark or antiquark line to the operator. As explained in 
PS, with the unimproved staggered action and unsmeared links, Y diagrams do not lead to mixing between different 
bilinears, and the result depends only on the "distance", A = X^^(<S' — F)^, of the bilinear. It is straightforward, 
though tedious, to check that these arguments generalize to the improved actions considered here. The result is 

1a = Ya-1 +Ia (A > 0) , (24) 

with Ya=o = 0, and 

Ia= [ BFsis^Pi'VY{A)+ [ 12BFsls2S^O^^VYiA) (25) 

Jk Jk 

1 tO\ C2 + C3 + Ci C2C3 + C2Ci + C3C4 

Vy{1) = 1 , Vy{2) = , Vy{3) = , Vy{A) = C2C3C4 . (26) 



The new functions are defined by 



P;:'ik) = D;:{k)D,{k) + J2 G.Akf , (27) 



and 



As,s,0;^^{k) ^ D;^G^,, + D,G,,,+ J2 (28) 



The single superscript 'W refiects the fact that the Naik term appears only at the quark-gluon vertex and not in the 

•)NN (-)N 



operator. Note that, unlike 0^1^ , Of^ is not symmetrical. 



C. Tadpole diagrams 

Here we include tadpole diagrams both on the external quark and antiquark propagators (i.e. self-energy contribu- 
tions), and those coming from the bilinear. In the latter case we include all diagrams in which the two gluons both 
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come from the bilinear, irrespective of whether they emanate from the same smeared Unk. Thiis, for the example of 
a distance-2 bilinear, which involves an average of a sum of products of two links, the gluon can couple between the 
links (as well as going from each link back to itself). 

It is convenient to divide the contribution into two parts, 

Ti = T^ + Tl, (29) 

with the former coming from ghion loops beginning and ending on the same smeared link, and the latter involving 
gluons propagating between smeared links. In both cases the result depends only on the distance A of the bilinear. 
No simple general formula covers all choices of links and action, so we quote the results in turn. 

1. For unimproved staggered fermions, the result is (PS) 

TX = (A-l)(^7r2-^i?/2) , (30) 



where the factor of tt^ comes from tadpole improvement using the fourth-root of the plaquette. If one uses the 
trace of the Landau gauge link, then tt^ is replaced by Jj. 3B/8. 

2. For the Fat-7 links (without the Naik term), the result is the same as for unimproved staggered fermions, 
eq. (30), due to cancellations. 

3. For O(a^) improved links (but without the Naik term), we find 



Tl = (A - 1) 



B/2 



2 / BciS2 



(31) 



4. For HYP links, we find 

n = {A-l) [ {-B/2)P, , (32) 

Jk 



where contains no Naik vertices: 



P^{k)=Df,{k)^ + J2G^Akf- (33) 



We emphasize that, at this stage, there is no tadpole improvement factor for the HYP links (although a related 

mean- field improvement will be introduced in sec. V). It is also noteworthy that this result would apply for 
both the Fat-7 and O(a^) improved links were one to also include projection back into SU(3) in those cases [24]. 



5. Finally, for the Asqtad-like action we find 



= (A - 1) 



(34) 



in which the second contribution is due to the Naik term. 

Now we turn to the "off-diagonal" tadpoles. These only arise from the bilinears, and not from the self-energy 
contributions, and are only present for operators with A > 2. Since they are off-diagonal they are not affected by 
SU{3) projection, and so take a common form for all actions and operators: 

Ti = 4 / Bs140i2Vt{A) (35) 

Jk 

Vt{2) = 1, Vt{3) = 2 + C3, Vt{3) = 3 + 2c3 + C3C4 . (36) 
Here O12 does not contain Naik contributions, even for the Asqtad-like action: 

4s^SpOw(^) = D^G^,p + DpGp,^+ J2 G-,>^G^,p. (37) 
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D. Self-energy diagrams 

If the Naik term is not included in the action, the flavor-singlet vector current is conserved, and it follows from the 
corresponding Ward Identity that its matching factor vanishes. Thus it must be that, in cases 1-4, 

Z = -Xu-Y^, i= (7^ <S) 1) . (38) 

Here we have used the result that Ti = for cases 1-4. Wc have checked eq. (38) analytically and mimcrically. 

We cannot use this relation for the Asqtad-like action, since the hypercube vector operator is not the conserved 
current (as it does not contain a Naik-like contribution). A direct calculation is needed, and we find: 

Z = j j BFIz , (39) 

Iz = ci(l + d^sj/2){l - 2{s^rF) {c?,Pr - Sc^Pi^^) 

- sisf (Pi + dN{2sl/i - cl)Di - dNslclci/Q) 

- 12sisf ci(l + dNsl/2)s2S^FO^^^ 

- 3s^S2S^ (4O12 - dNCiclGi,2 + dNslG2,i/Qj . (40) 

The first term in Z is the standard integral used to subtract the divergent piece. We have inserted rfjv in appropriate 
places so that this result is valid for all the actions we consider. 



V. FURTHER MEAN-FIELD IMPROVEMENT OF OPERATORS 



It is possible to apply another level of tadpole, or, more accurately, mean-field, improvement to the operators and 
actions which involve smeared links. Actually, for the HYP smeared links, this is the first level of tadpole improvement. 
The fluctuations in the smeared links are reduced compared to those of the original links but are still present. The 
residual fluctuations can be estimated and partially removed by defining a smeared mean-link by analogy with the 
definition of the original mean-link [18]: 

(ti^^)4 = (Smeared - Plaquette) . (41) 

Here the "smeared-plaquette" means the plaquette built out of smeared links. The operators are then mean- field 
improved by multiplying them by 

(no^^)i-^, (42) 

where A is the number of links in the bilinear. The argument leading to this factor is identical to that used in PS when 
tadpole improving staggered operators, and we do not repeat it here. This procedure should be simple to implement 
in practice. 

Wc have calculated the effect of such a mean-field improvement for Fat-7, O(a^) improved, and HYP links. We 
have not applied it to the case of the Asqtad-like action, because it is not entirely clear to us how to incorporate the 
Naik term. We find that one must add to the tadpole contribution the following: 



Tl = (A - 1) 



(43) 



where T^' and Pi are given in the previous section. The quantity in square parentheses evaluates to —0.9125 for 
Fat-7 links, -4.0634 for 0{a^) improved links, 0.5782 for HYP links with ai_3 = 0.75, 0.6, 0.3, and to 1.0538 for HYP 
links with the "Fat-7 choice" ai_3 = 7/8,4/7, 1/4. These values are substantially smaller than the analogous factor 
in the tadpole improvement of the unsmeared links, namely tt^. They are, nevertheless, significant, as we see in the 
next section. 

It is noteworthy that the Fat-7 and O(a^) improved smeared links receive a mean-field correction of opposite sign 

to that of both the HYP links. This indicates that the fluctuations in the former case have been "overcompensated" 
by smearing, and suggests that this higher level of mean-fleld improvement is likely to be more significant for the 
HYP smeared links. 

Finally, we note that after this higher level of mean-field improvement, the results for Fat-7 and HYP links with 
ai_3 = 7/8,4/7, 1/4 are identical. The equality of the tadpole contributions can be seen by combining eqs. (30), (32) 
and (43); that of other contributions follows from the fact that the single-gluon vertex is the same in both cases. 
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VI. NUMERICAL RESULTS AND DISCUSSION 



We present numerical results for the matching coefficients in Tables I-III. As explained in Ref. [17], the corrections 
are unchanged if the operators arc multiplied by (75 (8" ^5), due to the conserved axial symmetry. Thus we show results 
for only half the operators. Recall that we have chosen the NDR scheme {MS with an anticommuting 75) and set 
H= 1/a. We expect this to be a reasonable choice for the matching scale, but, in any case, the dependence of the Cij 
on /i is weak, as can be seen from eq. (IS).** 

The most striking result from the tables is the significant reduction in the size of one-loop corrections for all of 
the choices of smeared links. This is true also for the off-diagonal matching constants, although here the corrections 
were small to start with. We also see that the mean-field improvement of Section V leads to a significant further 
reduction in the corrections for HYP smearing, although the corrections increase somewhat for Fat-7 and 0{a^) 
improved smearing. 

To compare the different alternatives for improvement we quote, in Table IV, the range of variation of the diagonal 
coefficients Cu, both for a given spin (varying the fiavor), and for all spins and flavors. The range for a given spin is 
independent of the rcnormalization scale /i (since di is the same for all flavors) , and thus is a good measure of the size of 
lattice contributions to matching factors. The range for all spins and flavors does depend on n, but only rather weakly. 
The table shows that, of the alternatives we have compared, Fat-7 links, with or without mean-field improvement, 
and mean-field improved HYP or Fat-7 links lead to the smallest range of corrections. A similar conclusion holds if 
we consider the maximum magnitude of the corrections rather than the spread. 



Operator 


(a) 


(b) 


(c) 


(d) 


(e) 


(/) 


(1® 1) 


-29.3551 


1.8696 


-4.3917 


-2.1750 


-0.5939 


-0.0966 


(1®^m) 


-8.6416 


2.4633 


-2.5643 


-0.3301 


1.8394 


2.4633 




0.5657 


2.8990 


-2.8420 


-0.7999 


4.0139 


4.8653 


(1®^m5) 


5.2378 


3.3351 


-4.0469 


-2.1427 


6.0380 


7.2676 




8.7493 


3.7704 


-5.5793 


-3.7774 


7.9837 


9.6693 


(7m ® 1) 


0.0000 


0.0000 


0.0000 


1.4155 


0.0000 


0.0000 


(Tm ® Cm) 


-4.9092 


0.7869 


2.9240 


4.2755 


-0.9457 


-1.1794 




0.1721 


-0.1201 


-2.9799 


-1.5110 


1.3090 


1.8461 


(7m ® Cm-') 


-3.3948 


0.3636 


-0.0621 


1.4290 


0.2617 


0.3636 


(7m ® C^m) 


2.5040 


-0.1930 


-5.4907 


-4.0295 


2.7140 


3.7396 


(7m ® Ci-s) 


0.1902 


0.1367 


-2.5010 


-1.0264 


1.6009 


2.1030 


(7m ® Cms) 


4.8930 


-0.2147 


-7.9437 


-6.4957 


4.1592 


5.6841 


(7m ® 6) 


2.7709 


0.0369 


-5.0332 


-3.5799 


2.9898 


3.9694 


(7m- ® 1) 


1.5969 


0.3741 


-1.3115 


-0.0393 


1.9442 


2.3404 


(7m-' ®Cm) 


0.8194 


0.8758 


2.1260 


3.3442 


0.9819 


0.8758 


(7m-' ® Cp) 


3.0150 


0.0410 


-4.4862 


-3.1752 


3.0313 


3.9735 


(7m'- ® Cm-) 


4.5728 


1.7594 


6.6960 


7.7590 


0.2703 


-0.2069 


(7m-' ® Cmp) 


1.2809 


0.3800 


-1.4041 


-0.1178 


1.9309 


2.3463 


(7m-' ®Cp<^) 


4.9409 


-0.2098 


-7.3985 


-6.0684 


4.2177 


5.6890 



TABLE L Diagonal part of the one-loop matching constants, ca, using the NDR scheme with /iNDR = 1/a in the continuum. The 
components fi, u, p and a arc all different. Results are given for six choices of action and operators: (a) unimproved; (6) Fat-7 links; (c) 
fully O(a^) improved links; (d) fully 0{a^) improved links and Naik term (Asqtad-like action); (e) HYP links with the smearing coefficients 
from Ref. [16] , ai-3 = 0.75, 0.6, 0.3; (/) HYP links with tree-level improvement coefficients, ais = 7/8, 4/7, 1/4. The error in the results 
is no larger than 0.0001. 



Approximate methods of calculating the optimal matching scale, q* , do not obviously generalize to the case of operators with 
non- vanishing anomalous dimensions. We return to this issue elsewhere [24]. 
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Name 


Operator-i 


Operator-j 


(a) 


(6) 


(c) 


(d) 


(e) 


CVVM 


(7m ®^-) 


(7m ® Cm) 


3.0412 


0.3508 


1.4104 


1.2976 


0.4203 


CVAM 


(7m ® Cms) 


(7m ® 6^5) 


-0.6463 


-0.2565 


-0.6192 


-0.5481 


-0.3000 


CVTM 


(7m ® 




-1.4861 


-0.2797 


-0.9241 


-0.8211 


-0.3512 


(-'7'_1..A/ 






-O.fiTG:! 




-().2()r)o 


-0.175:! 


-0.0202 



TABLE II. Non vanishing ofT-diagonal one-loop matciiing constants, Cij. The components fi, v and p are all different, but otherwise 
can take any values. Results are given for same choices of action and operators as in Table (I), except that HYP links with tree- level 
coefficients give identical results to Fat-7 links ajid thus are not shown. The error in the results is no larger than 0.0001. 



Operator 






(&',/') 






(c') 






(e') 


(1® 1) 






0.9571 






-8.4551 






-0.0156 


(1 ® ?m) 






A coo 

2.4DOO 






v^a AO 






1.8394 


(1®Cm-) 






3.8115 






1.2214 






3.4357 


U OS) ?M5j 






5.1600 






4.0799 






4.8815 


(I® 6) 






6.5079 






6.6110 






6.2490 


(7m ® 1) 






0.0000 






0.0000 






0.0000 


(7m ® Cm) 






-0.1255 






-1.1394 






-0.3675 


(7m ^C-) 






0.7924 






1.0835 






0.7308 


(7m ® Cm^) 






0.3636 






-0.0621 






0.2617 


(7m ® C^p) 






1.6320 






2.6361 






1.5576 


(7m ® Ci-s) 






1.0492 






1.5624 






1.0227 


(7m ® Cms) 






2.5227 






4.2466 






2.4245 


(7m ® Cs) 






1.8619 






3.0936 






1.8334 


(7m- ® 1) 






1.2866 






2.7519 






1.3660 


(7m-' ®Cm) 






0.8758 






2.1260 






0.9819 


(7m'- ® Cp) 






1.8659 






3.6407 






1.8748 


(7m'^ ® Cm--) 






0.8469 






2.6325 






0.8486 


(7m-' ® Cmp) 






1.2925 






2.6593 






1.3527 








2.5276 






4.7918 






2.4830 


TABLE III 


. Results for 


cu after the 


mean-field improvement discussed in sec. V. Notation as in 


Table (I). 


Results are given for (6', /') 


Fat-7 links and HYP links with tree-level improvement coefficients, ai- 


-3 = 


7/8,4/7, 1/4; (c') fully 0{a?) improved Unks; (e') HYP links 


with smearing coefficients from Ref. [16], ai_3 = 0.75, 0.6, 0.3. 












Spin 


(a) 


(&) 


(c) 


(d) 


(e) 


(/) 


(&',/') 


(c') 


(e') 


S/P 


38.1 


1.9 


3.0 


3.4 


8.6 


9.8 


5.5 


15.1 


6.3 


V/A 


9.8 


1.0 


10.9 


10.8 


5.1 


6.9 


2.6 


5.4 


2.8 


T 


4.1 


2.0 


14.1 


13.8 


3.9 


5.9 


1.7 


2.7 


1.6 


All 


38.1 


4.0 


14.6 


14.3 


8.9 


10.8 


6.6 


15.1 


6.6 



TABLE IV. Spread of values for diagonal corrections, ca, both for a given spin, and between all operators. Notation as in Tables (I) 
and (III). 
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What values of ca give rise to "small enough" corrections in present simulations? Taking l/o, = 2 GeV as a typical 
lattice spacing, and using Q:p^(2GeV) ~ 0.19, we find Cpaj^/{A7r) « 0.02. Thus a matching coefficient c = 5 
corresponds to about a 10% correction at this lattice spacing. This is the size of corrections we are aiming for, and 
we see that tadpole improved HYP fermions lead to corrections of about this size. 

Finally, it is interesting to compare to the size of one-loop corrections for bilinears obtained with other fermion 
actions. For unimproved Wilson fermions one finds, after tadpole improvement (picking for definiteness the tadpole 
improvement scheme of Ref. [25]), Cj = —0.1, —9.7, —7.8, —2.9, —4.3 for i = S,P, V, A, T, using the same renormaliza- 
tion scheme and scale for the continuum operator as in the tables. We have not been able to find the corresponding 
results for improved Wilson fermions incorporating tadpole improvement, but it is clear from Table 3 of Ref. [26] 
that one-loop matching factors are of similar size as for unimproved Wilson fermions. For domain-wall fermions, 
tadpole-improved results are given in Ref. [27]: Cj = -11.2,-5.3,-2.0 for i = S/P,V/A,T (setting the domain-wall 
mass M = 1.7). We conclude that the size of corrections with improved staggered fermions is comparable to, or 
smaller than, that for other fermions. This provides further impetus to pursue calculations with improved staggered 
fermions. 
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APPENDIX A: NOTATION AND TECHNICAL DETAILS 



We use the hypercube fields and bilinears introduced in Ref. [28] . The lattice is divided into 2^ hypercubes labeled 

by a vector y, with all components even. Points within a hypercube arc labeled by a hypercube vector (C^, -D^ in the 
following), with all components or 1. The lattice bilinears we use are specified by "spin" and "flavor" hypercube 
vectors and in the following way in terms of the staggered fields x and x- 

Oi^s^i,) (y) = 4 E + cKis^Cf^dx' {y + D). (ai) 

CD 



Here the sums run over all positions in the hypercube, and 

1, 



{is ® ^f)cd = [rclslDYFl , (A2) 
with 

75 = ltli'l3'lt (A3) 

composed of hermitian Euclidean gamma matrices. It follows that the x and x fields are separated by a fixed number 
of links which is given by the "distance" A = X]^(S'^ — ^m)^- continuum limit, this lattice bilinear has the 

same spin, flavor and normalization as the continuum bilinear 

Ol^lM=QlalfeF'Qh, (A4) 

where 6 is a four-flavor quark field, with spinor index /3 and fiavor index b, both running from 1 to 4, and the 
^F = "fp form a convenient basis for the flavor matrices. 

The superscript on x*^, Q*^, etc. indicates an additional flavor index, corresponding to the different continuum 
flavors (u, d, s, etc.). We consider here only continuum flavor non-singlet operators, so that we do not have to 
calculate diagrams in which the fermion fields in the bilinear are contracted together. 

Because the lattice bilinears are spread over a hypercube, the phase factors due the external quark and antiquark 
momenta depend upon the hypercube vectors C and D in eq. (Al). To disentangle these phases Daniel and Sheard 
introduced the following definition [29]: 



ne"=^/ME^M(fc)(-)^-^ (A5) 



M 



with 



EM{k) =l[l (e-''=''/2 + (-)^^e''=''/2) . (A6) 

Here we have introduced the conjugate hypercube vectors 

=2 E (A7) 

Previous results for X diagrams in Ref. [29] and PS were expressed as sums over the additional hypercube vector M, 
which were then evaluated numerically. We have been able to perform this sum analytically, which simplifies the final 
expressions. The two results we need in this paper are: 

V,{k) = Y,EM{k)EM{-k){-)^-(^+'^^ =J{cos[k^{S-F)^] , (A8) 

M n 

where i labels the spin and flavor, and 

Y,EM{k)EM+^.+.{-k){-)^-(^+'^^ = -\{S~F)^\\{S-FUs^s^cos[kp{S-F)p]cos%{S-F)„], (A9) 

M 
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in which all indices arc different. These results can be obtained by combining PS eqs. (24) and (27). 
The explicit forms for lattice integrals are abbreviated using the following notations: 

f dk 

J = 167r^ Yi J ^ ' V = sin(fc^/2) , = sin(fc^) , = cos(A;^/2) , = cos(A;^) . 
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